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Current experimental results suggest that some organic quasi-one-dimensional superconductors 
exhibit triplet pairing symmetry. Thus, we discuss several potential triplet order parameters for the 
superconducting state of these systems within the functional integral formulation. We compare weak 
spin-orbit coupling fxyz, Px, Py and pz symmetries via several thermodynamic quantities. For each 
symmetry, we analyse the temperature dependences of the order parameter, condensation energy, 
specific heat, and superfluid density tensor. 
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Recent NMR experiments in the organic conductor 
(Bechgaard salt) (TMTSF)2PF6^'^ in combination with 
earlier upper critical field (-ffc2 ) measurements'^ indicate 
that this system may be an unconventional triplet su- 
perconductor at low temperature (T). Lee et. almi^ 
observed that there is no ''^Se Knight shift for fields 
H II b' (P « 6 kbar) and H || a (P « 7 kbar). This indi- 
cates that the spin susceptibilities in the superconduct- 
ing state are Xa ~ Xn, and Xf ~ Xn, where xn is the 
normal state susceptibility. Their measurements suggest 
that the d-vector order parameter for triplet supercon- 
ductivity is pointing along the c* axis. Unfortunately, 
the spin susceptibility in the superconducting state for 
fields along the c* axis could not be measured because 
the upper critical field {H^^ ) was exceeded before a sig- 
nal could be detected. Since the low T dependence of 
Xc* (T) is not accessible experimentally, it may not be 
possible to use this technique to decipher the node struc- 
ture of the triplet order parameter. The symmetry of the 
order parameter in the sister compound (TMTSF)2C104 
was preliminarily explored in the thermal conductivity 
experiments of Belin and Behnia* (BB). They indicated 
that their data was inconsistent with the existence of gap 
nodes at the Fermi surface as suggested by Takigawa et. 
al^. These results combined with upper critical field 
measurementsSii seem to suggest a fully gapped triplet 
state in (TMTSF)2C104. However, detailed ^^Se NMR 
experiments seem to be lacking for (TMTSF)2C104. 
These experiments were inspired by early suggestions of 
triplet superconductivity in the Bechgaard salta^iSiifliii, 
but following these new experimental developments the- 
oretical efforts intensified -^^i^^i^^^i^^i^^. Lebed, Machida, 
and Ozaki (LMO)^^ proposed a "p-wave" triplet order 
parameter for (TMTSF)2PF6, where the d-vector had a 
strong component along the a direction, thus producing 
a strongly anisotropic spin susceptibility with Xa ^ Xn 
and Xb' ~ Xn- A fully gapped singlet "d-wave" order 
parameter for (TMTSF)2C104 was proposed by Shima- 
hara^'^, while gapless triplet "f-wave" superconductiv- 
ity for (TMTSF)2PF6 was proposed by Kuroki, Arita, 
and Aoki (KAA)i4. Duncan, Vaccarella and Sa de Melo 
(DVS)i^*i& performed a detailed group theoretical anal- 
ysis and suggested that a weak spin-orbit fully gapped 
triplet "pa;-wave" order parameter, where Xa ~ Xn and 



Xb' ~ Xn^^^^, would be a good candidate for supercon- 
ductivity in Bechgaard salts. 

In this paper, we consider only triplet states correspond- 
ing to weak spin-orbit coupling, since in (TMTSF)2X 
the heaviest element is ^^Se. Within the orthorhombic 
{D2h) group, this limits the number of possibilities for 
the order parameterii^ii^ to four symmetries, Px, Py, Pz, 
fxyz- For each one of these symmetries, we calculate the 
temperature dependences of the order parameter, con- 
densation energy, entropy, specific heat, and superfluid 
density tensor. 

We study single band quasi-one-dimensional systems in 
an orthorhombic lattice with dispersion relation 

Ck = - \tx\cos{kxax) - \ty\cos{kyay) - \tz\cos{kzaz), 

where \tx\ S> |tj,| S> |tz|- Furthermore, a^, ay and 
in our notation correspond to the unit cell lengths along 
the crystallographic directions a, b', and c* respectively. 
We work with the Hamiltonian H = Hkm + Hint , where 
the kinetic energy part is Hkin = J2k,a £.ktpk,a^Ka, with 
= Ek — M and the interaction part is 

H^nt = ^ E E VcffysiK ^')bUK q)^5(k', q) (1) 

kk'q a/S-yS 

with blp{k,q) = V'ik-Hq/2,aV'k+q/2,/3' ^^^"^^ Z^' 7 and 
5 are spin indices and k, k' and q represent linear mo- 
menta. We use units where h = kg = 1. 
In the case of weak spin-orbit coupling and triplet pair- 
ing, the model interaction tensor can be chosen to be 

Fa^^5(k,k') = Vt/ir(k,k')0r(k)0Mk')r (2) 

where Ta/^^s = Vq/j • with v^^ = (iaay)ap. Vr 

is a prefactor with dimensions of energy which char- 
acterizes a given symmetry. Furthermore, the term 
/ir(k, k')0r(k)(/)J(k') contains the momentum and sym- 
metry dependence of the interaction of the irreducible 
representation F with basis function 0r(k) and 0p(k') 
representative of the orthorhombic group {D2h)- 
We use the functional integration method to write down 
the partition function Z = J I? [?/;^, -0] exp [5] where 
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S = JdT [Ek,a<aM(-5-)V'k,o(r) -ff(Vt,^)J . We 

treat Hi„t in the zero center of mass momentum 
(q = 0) pairing approximation. Furthermore, since 
we are mostly interested in symmetry aspects we 
take /ir(k, k') — 1. We introduce the Gaussian in- 
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exp[-(5G], where Qg 



I dr Ekk' ^1 (k, (k', r) /Vp. The shift transforma- 
tion I?|(k,T) ^ I?J(k,r) + yrVa/3,»0r(k)V^Lk,„V'k,^ elim- 
inates Hint and integration over the fermionic degrees of 
freedom results in the effective action 

^•cff = -QG-Tr In [M/2], (3) 

where M is a 2 x 2 block diagonal matrix of the form 

[dr + Ck] Sap Y.k' 'r)v/3a.i(/'r(k) 

Ek' Vlik', r)vl^ >f (k) [-a. - Ck] S^p 



Thus, Z = Jd vl^Vi exp[S'off]. At the saddle point 

approximation ^^^(k, r) is taken to be r independent, 
and I?i(k') = %Ar. The order parameter equation 

is obtained from the stationary condition 5S^^ / Svl = 
leading to 



J2 Vr\^rO<.)\hanhiEk/2T)/2E^, 



(4) 
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where £^k = a/^^ + l^rp|0r(k)p. The number equation 
is obtained from N = —dflo/d^, where fla = —TS^,^ is 
the saddle point thermodynamic potential, and results 
in 



N 



k 



(5) 



where rik = [1 — ^ktanh(iJk/2r)/£^k] is the momentum 
distribution. These two equations must be solved self- 
consistently, and quite generally they are correct even 
in the strong coupling (or low density) regime provided 
that T <^ Tc- Corrections to and N can be ob- 
tained by considering Gaussian fluctuations. Writing 
X',(k,r) = X',(k) + (5X',(k,r) and expanding rrln[M] 
to quadratic order in SVi (k, t) results in the effective 

action 5off = sf^ - Tr [M-^U] ^ /2, where U is a 2 x 2 
block matrix that contains only off-diagonal elements, 
with Ui2 = X^k' '5^«(k','r)v/3a,j'/'r(k) and U21 = UI2. 
These corrections are important for small carrier den- 
sity (or Tc ~ Ep), however, in the BCS limit of high 
carrier density (or Tc <C Ep), discussed here, these cor- 
rections are negligible for N, but are important for Q 
only when T < Tc. For Bechgaard salts, Tc w 1.5 K 
and Ep « 3,071 K (with respect to the bottom of 
the band). The vector T>i is related to the standard 
d-vector via the relation (ii(k) — ^j^, I?i(k')^r(k). In 
the D2h point group all representations are one dimen- 
sional and non-degenerateiS*i&, which means that the d- 
vector in momentum space for unitary triplet states in 
the weak spin-orbit coupling limit is characterized by one 



FIG. 1: Temperature dependence for (a) Ar and (b) /ir in 
thermal units (K). Notice that /ir is nearly T independent. 



of the four states: (1) ^^i„(a), with d(k) = TjAf^^^XYZ 
C'f.yz" State); (2) ^Bi^{a), with d(k) = r^Ap T ("p," 
state); (3) ^B2u{a), with d(k) = fjAp^Y (%" state); 
(4) 3^3^,(0), with d(k) = vAp^X ("p^" state). Since, 
the Fermi surface touches the Brillouin zone boundaries 
the functions X, Y, and Z need to be periodic and 
can be chosen to be X = siii{kxax), Y = sin{kyay), 
and Z = sin(fc2az). The unit vector 77 defines the di- 
rection of d(k). The T dependence of Ap and fir for 
fxyz, Px, Py and Pz symmetries are shown in Fig. ^ 
The parameters used are \tx\ = 5800 K, \ty\ = 1226 K 
and \tz\ = 48 K, and N/Nmax = 1/4 (quarter-filling). 
The minimum value of the dispersion e(k) is tmin = 
—7,074 K, while the maximum is emax = +7,074 K. 
In order to make direct comparison between different 
symmetries we choose Tc = 1.5 K for all symmetries. 
This requirement forces the interaction strength pref- 
actors to be Vf^^^ « -12,577 K; Vp, « -2,910 K; 
Vp^ « -3,025 K;Vp. ~ -3,208 K. Notice in Fig. Uti 
that Af^^^ > Ap^ > Ap^ > Ap^ for all T < Tc, and 
that fir is largely independent of T, however, the condi- 
tion fLp^ > fip^ > fif^^^ > fip^ applies for all T < Tc. 
The relative stability of these phases can be studied 
via the condensation energy ATr = ^r{T) — J-n{T), 
where Tr{T) and Tn{T) are the Helmholtz free ener- 
gies of the superconducting state with symmetry F and 
of the normal state, respectively. ATt is calculated us- 
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FIG. 2: Temperature dependence of the condensation energy 
AjTr in thermal units {K). 



FIG. 3: Temperature dependence of specific heat Cr in units 
of ks- Inset shows low T behavior. 



ing = + /iiV. Fig. |21 shows AjFp for the fxyz, Px, 
Py and Pz symmetries. Notice that AJ-'r are very small, 
are expressed in thermal units, and obey the relations 



> \AJ^f 



for T <Tr. The 



Px symmetry which has a full gap, has also the largest 
(negative) ATy. The values of ATp, and ATp^ are very 
close reflecting that the and Py symmetries have lines 
of nodes at Z = and Y — 0, respectively. Lastly, the 
fxyz symmetry has double zeros ai Z — and F = and 
has line nodes at Z = or F = 0, which are costly in 
condensation energy. This confirms the general expec- 
tation that a fully gapped (nodeless) superconducting 
phase {px) is more likely to win over competing phases 
which have nodes {py, Pz and fxyz)- 
Since fig is only a function of |d(k)p (recall that 
|d(k)|2 = |Arp|(/)r(k)|2) the vector nature of d(k) 
does not appear explicitly in thermodynamic properties. 
However, the node structure of triplet order parameters 
can be probed by a specific heat measurement. The spe- 
cific heat Cr = -Td^J'r/dT^ is 



2 ^^^ Tdmi 



where P{E^^) = f {Ei,) [I ~ f {Ey,)] . The results for 
Cr can be seen in Fig. and analysed as follows. 
The specific heat jump at Tc = 1.5 K is character- 
ized by the parameter Or — Csx{Tc)/Cn{Tc) — 1, 
which takes values 6*/^^^ = 0.5596, 9p^ — 0.9733, 
Op^ = 0.9358, and Op^ = 1.2001. The corresponding 
singlet s-wave value is Og = 1.4604. An analysis of Cr 
at low T is also important in distinguishing possible 
weak spin-orbit triplet phases. Using the fact that 
dfi/dT « 0, d\d{k}\^/dT « at low T, then Cr « 
2T dww'^sech'^{w)N{2Tw), where w = lj/2T and 
N{uj) — 2 J2it: S{uj — E]i) is the auxiliary density of states. 



This leads to 
where the 



C 



iT/Tf^J^log{K,Af^^jT), 



temperature scale 
2{Af^^Jxy^ D/ttk 7(3), with D = 
7(to) = dw'w"^sech^{w) 

Kr, 



K 
2/ 



T7^ 

Jxyz 

16/[7r|sin(a;)|], 
and 



— exp [Jp°° dww™sech^(u;)log [4w^-^k] /^{m)] 



This T dependence of Cf^^, is a reflexion of the the 
behavior of N{uj) — a(tj/Ao)log(/3Ao/w) at small 
Lo, where a and /3 are constants. The presence of 
the logarithmic term results from the large number 
of low-energy states surrounding the double zeros 
of the order parameter at the relatively flat Fermi 
surface. The Pz and Py symmetries, however, have 
very different behaviors at low T, given that both 
order parameters have lines of nodes at the Fermi 
surface, but no double zeros. In this case, the specific 



Py 



(T/Tp 



heats behave as Cp, — (T/Tp^)'^, and Cj 

= [Apjixl] ^£'7(3), for {pj = Pz or py) 



where T„ ^ 



Ap 
with 



Lastly, Cp^ = (TpjT)^/\xpi-LUo/T), where loq = 
\ty\ + \tz\- fi) /\tx\f and TpY' - ^^7(1/2), 



1 - 
A 



2iu;oy/V[7T^\sinix)\Al{\ty\ + \t^\-f,)], 



B = [{tl - AlJ/{tytz)]^^'^ . Notice that Cp^ is ex- 
ponentially supressed at low T due to the presence 
of a full gap in the excitation spectrum, however the 
exponential prefactor behaves as T"^/^, which differs 
from the T~^/^ behavior of the prefactor in the singlet 
s-wave case. This difference arises because N{uj) in the 
Px case has a square root dependence in the vicinity of 
the gap edge, Np^ ^ ^uj — wo^p^ , while in the s-wave 
case N{uj) has a square root singularity near the gap 
edge A^s ~ 1/ — wo,s- The Gaussian correction fie to 
fio leads to collective modes at low T with anisotropic 
dispersions uj = cr(0, 0)|q| for all symmetries. These 
modes give a symmetry independent contribution 
to Cr, but the prefactor is symmetry dependent. This 
collective mode contribution is characteristic of neutral 
superfluids. However, in real charged superconductors 
that may become plasmonized, and thus gapped. 
The temperature dependence of the superfluid density 
tensor pij (T) can also be used to distinguish different 
weak spin-orbit triplet phases. This tensor is directly 
associated with phase twists of the U{1) phase of the 
d-vector. Take d(k) d(k)exp [i0(k)] and expand Scs 
in powers of (/'(k) about the saddle point with 0(k) = 0. 
The resulting action AS" = S'off (0) - Scs{(l) 0) is AS" = 
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-V/2Y.(l){k)(f,{-k)hkjp,j, with 

Pzj{T) = X! ['^k^i^j^k - ikSi^k^jCk] , (6) 

k 

where is the momentum distribution, and 

yk = {2T)-'^sech^ {Ek/2T) is the Yoshida distri- 
bution. In Fig. 0] we show the T dependence of pij . 
In the case of the I?2/i group only diagonal compo- 
nents Pa exist, but they are highly anisotropic due 
to the quasi-one-dimensionality of ^k- The low T 
behavior of Apu = [pii{T) / pii{0) — 1] is shown in the 
inset of Fig. |3| At low T, the main contributions 
to Apii come from the term containing Ik- For the 
symmetry, Ap,, = ~iT/Tiy/\xp{-LUo/T), 

-(T/Tjfj3/2exp(-tJo/r), where T^^ are character- 
istic temperatures. Notice the exponential behavior due 
to presence of a full gap. Furthermore, notice that the 
T dependence of the prefactor of Ap^x is different from 
Apyy and Ap^z due to the highly anisotropic Fermi 
surface of these systems, i.e., the velocity = dx£,i<i 
does not vanish anywhere at the Fermi surface. For 
the py symmetry, Ap^x = -T/T^^, Apyy - -{T/T^ f, 
and Apzz = —T/T^ . Notice the power law for 
Apyy, which results from the simultaneous contribu- 



tion from the lines of nodes of dp^ and the zeros of 
Vy = i9j,^k- Similar behavior is found also for the 
Pz symmetry, Apxx = -{T/Tl), Apyy = -{T/Tyj, 
and Apzz — —{T/TpJ^. In the Pz case, however, 
the dependence appears in Apzz- Finally for the 
]xy. symmetry, Apxx = -(T/T- ^_)log(i^iA/^„ JT) 
Apyy = -T/Tl^^, and Apzz = -T/Tl^^. The loga- 
rithmic dependence on Apxx originates from the node 
structure of d^^ ^ and from a non- vanishing Vx at the 
Fermi surface. Notice that the logarithmic dependence 
is absent on Apyy and Apzz, because the zeros of Vy 
and Vz cancel it out. The T dependence of pii(T) can 
be measured via penetration depth experiments. 
In summary, we studied several properties of organic 
quasi-one-dimensional conductors (Bechgaard salts) 
which are strong candidates for triplet superconductiv- 
ity. We compared weak spin-orbit coupling symmetries 
fxyz, Px, Py and Pz via several thermodynamic quanti- 
ties. For each symmetry, we analysed the temperature 
dependences of the order parameter, condensation en- 
ergy, specific heat, and superfluid density tensor. We 
would hke to thank NSF (Grant No. DMR-9803111) 
for support. ^Present address: Department of Physics, 
Massachusetts Institute of Technology, Cambridge, MA 
02139-4307. 
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FIG. 4: Temperature dependence of Vopu/ai in thermal 
units {K). Vo (ttj) is the unit cell volume (length). Insets 
show the low T behavior of Apu = [pu (T) / pu (0) — 1] . 



